Heavy sfermions in SUSY analysis at LHC and ILC 



K. Desch\ J. Kalinowski^, G. Moortgat-Pick^^, K. RolbieckP, and W.J. Stirling^ 

^ Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 

Instytut Fizyki Teoretycznej, Uniwersytet Warszawski, PL-00681, Warsaw, Poland 
^ IPPP, Institute of Particle Physics Phenomenology, University of Durham, Durham DHl 3LE, UK 



Abstract. The physics potential of the Large Hadron Collider in combination with the planned 
International Linear Collider is discussed for a difficult region of supersymmetry that is charac- 
terized by scalar SUSY particles with masses around 2 TeV. Precision measurements of masses, 
cross sections and forward-backward asymmetries allow to determine the fundamental supersym- 
metric parameters even if only a small part of the spectrum is accessible. No assumptions on a 
specific SUSY-breaking mechanism are imposed. Mass contraints for the kinematically inaccessible 
particles can be derived. 

PACS. 14.80.Ly Supersymmetric partners of known particles - ll.30.Pb Supersymmetry 



1 Introduction 

Supersymmetry (SUSY) is one of the best-motivated 
candidates for physics beyond the Standard Model (SM) 
If experiments at future accelerators, the Large Hadron 
Collider (LHC) and the International Linear Collider 
(ILC), discover SUSY they will have as well to de- 
termine precisely the underlying SUSY-breaking sce- 
nario with as few theoretical prejudices as possible. 
Scenarios where the scalar SUSY particle sector is very 
heavy as required, for instance, in focus-point scenarios 
(FP) [1] are particularly challenging. In these scenar- 
ios the gaugino masses are kept relatively small while 
squark and slepton masses might be in the multi-TeV 
range. It is therefore of particular interest to verify 
whether the interplay of a combined LHC/ILC anal- 
ysis p] could unravel such models with very heavy 
sfermions. 

Methods to derive the SUSY parameters at col- 
lider experiments have been worked out, for instance 
in [3l|4]. In [5l[6ll7] the chargino and neutralino sec- 
tors have been exploited at the ILC to determine the 
MSSM parameters. However, in most cases only the 
production processes have been studied. Furthermore, 
it has been assumed that the masses of the virtual 
scalar particles are already known. In the case of heavy 
scalars such assumptions, however, cannot be applied 
and further observables have to be used to determine 
the underlying parameters. Studies have been made 
to exploit the whole production-and-decay process in 
the chargino/neutralino sector [8]. Exploiting spin ef- 
fects, it has been shown in [9] that, once the chargino 
parameters are known, useful indirect bounds for the 
mass of the heavy virtual particles could be derived 
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from forward-backward asymmetries of the final lep- 
ton Afb(^)- 

Here a FP-inspired scenario is discussed that is 
characterized by a ^ 2 TeV scalar particles sector [10] . 
In addition, the neutralino sector turns out to have 
very low production cross sections in e^e~ collisions, 
so that it might not be fully exploitable. Only the 
chargino pair production process has high rates and 
all information obtainable from this sector has to be 
used. The analysis is performed entirely at the EW 
scale, without any reference to the underlying SUSY- 
breaking mechanism. Masses, cross sections and spin- 
dependent forward-backward asymmetries are mea- 
sured at the LHC and at the ILC with y/s < 500 GeV. 
The potential of a multiparameter fit to determine the 
underlying parameters has been analysed. 

2 Case study at LHC and ILC 

We study chargino production with subsequent lep- 
tonic and hadronic decays 

e'+e+^xt+Xl, (1) 

xt^Xi+i^ + i^ and Xi+Qd + qu, (2) 

Xl^Xi+^'+i^ and Xi+qd + Qu, (3) 

where £ = e, ^, qu = u,c, — d,s. The produc- 
tion process contains contributions from 7- and Z*^- 
exchange in the s-channel and from zz-exchange in the 
i-channel. The decay processes are mediated by W^, 
£l, or by qdh, quh exchange. The masses and eigen- 
states of the neutralinos and charginos are determined 
by the fundamental SUSY parameters: the C/(l), SU{2) 
gaugino mass parameters Mi, M2, the Higgs mass pa- 
rameter ^ and the ratio of the vacuum expectation 
values of the two neutral Higgs fields, tan/3 = 
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The chosen MSSM scenario can be characterized 
via the following mSUGRA parameters, taken at the 
GUT scale except for tan/3: mi/2 — 144 GeV, mo = 
2 TeV, ^0 = GeV, tan/3 = 20, sgn(/i) = +1. How- 
ever, our analysis is performed completely within the 
general MSSM framework, without any reference to 
the underlying SUSY-breaking mechanism. The pa- 
rameters at the EW scale are given by 

Ml = 60 GeV, M2 = 121 GeV, M3 = 322 GeV 
Ai = 540 GcV, tan/3 = 20 (4) 

The derived masses of the SUSY particles are listed in 
Table [T] As can be seen, the charginos and neutralinos 
as well as the gluino are rather light, whereas the scalar 
SUSY particles have masses about 2 TeV. 

2.1 Expectations at the LHC 

All squarks in this scenario are kinematically accessi- 
ble at the LHC. The largest squark production cross 
section is for ti^2 - However, with stops decaying mainly 
to gt (with BR{ii,2 gt) ^ 66%), where background 
fr'om top production will be large, no new interesting 
channels are open in their decays. The other squarks 
decay mainly via gq, but since the squarks are very 
heavy, r ^ 2 TeV, precise mass reconstruction 
will be difficult. 

Since the gluino is rather light in this scenario, 
several gluino decay channels can be exploited. The 
largest branching ratio for the gluino decay in our sce- 
nario is a three-body decay into neutralinos, BR{g — > 
ji^hh) ^ 14%, followed by a subsequent three-body lep- 
tonic neutrahno decay BR{x2 Xi^^^ ); i — e,^ 
of about 6%, see Table [2j In this channel the dilep- 
ton edge will be clearly visible, since this process has 
low backgrounds [2]. The mass difference between the 
two light neutralino masses can be measured from the 
dilepton edge with an uncertainty of about [TT] : 

(5(m-^o — m^o) ~ 0.5 GeV. (5) 

The gluino mass can be reconstructed in a manner 
similar to the one proposed in [12], where the SPSla 
scenario is analysed. Although our scenario is different, 
in both we are systematics limited due to hadronic 
energy scale and a similar relative uncertainty of ^2% 
can be expected. 

2.2 Expectations at the ILC 

At the ffi-st stage of the ILC, ^/s < 500 GeV, only light 
charginos and neutralinos are kinematically accessible. 
However, in this scenario the neutralino sector is char- 
acterized by very low production cross sections, below 
1 fb, so that it might not be fully exploitable. The low 
cross sections are due to the mixing character of the 
neutrahnos, see pUj . 

Only the chargino pair production process has high 
rates. We constrain our analysis to the ffist stage of the 



ILC with y/s < 500 GeV and study only the xfxi 
production. 

The chargino mass can be measured at = 350 
and 500 GeV in the continuum, with an error of about 
0.5 GeV [Hill]. This can serve to optimize the ILC 
scan at the threshold which, because of the steepness 
of the s-wave excitation curve in XiXi production, 
can be used to determine the light chargino mass very 
precisely, to about [T4] : 

m^± = 117.1 ±0.1 GeV. (6) 

The mass of the lightest neutralino m^o can be de- 
rived via the decays of the light chargino, cither from 
the energy distribution of the lepton £^ {BR{xi 
X^l~i>i) ^ 11%, see Table [2]) or from the invariant 
mass distribution of the two jets in hadronic decays 
{BR{xi XiQdQu) 33%, see Table [2]). We take [E] 

m^o = 59.2 ± 0.2 GeV. (7) 

Together with the information from the LHC, Eq. ([5]), 
a mass uncertainty for the second lightest neutralino 
of about 

m^o = 117.1 ±0.5 GeV (8) 

can be assumed. 

We identify the chargino pair production process, 
e~^e~ XiXi 7 i'^ the fully leptonic (^^!^Xi^~^Xi) ^^'^ 
semileptonic {£i^XiQQ'Xi) final states (where i = e,fi). 
We estimate an overall selection efficiency of 50%. For 
both final states, W^W" production is expected to 
be the dominant SM background. For the semileptonic 
(sic) final state, this background can be efficiently re- 
duced from the reconstruction of the hadronic invari- 
ant mass. 

In Table [3l we list cross sections multiplied by the 
branching fraction B^ic = 2 x BR{xi XiQdQu) x 

Bi?(xr ^ Xi^'l^) + [BR{xi Xi^'i^)? ~ 0.34 (first 
two families) including an e^ic = 50% selection effi- 
ciency. The error includes the statistical uncertainty 
based on £ = 200 fb^^ in each polarization configura- 
tion, (Pe-,^e+) = (-90%, +60%) and (+90%, --60%), 
and a relative uncertainty in the polarization of AP^± / Pe± = 
0.5% [Hj. 

The statistical error on the forward-backward asym- 
metry ApB, based on binomial distribution, is given by 

Z\(Afb)^*"' = 2^e{l-e)/N, (9) 

where ctf.b are the acceptance-corrected cross sections, 
e = (Tf/((Tf + cb) and N denotes the number of se- 
lected events. In Table [3] the asymmetries are listed 
only for the {Pe-,Pe+) = (-90%, +60%) case, since 
the cross sections for the opposite polarization are very 
small and the statistical errors become very large. Con- 
sequently we do not include them in the following anal- 
ysis. 

3 Parameter determination 

We determine the underlying SUSY parameters in sev- 
eral steps: 



Desch, Kalinowski, Moortgat-Pick, Rolbiecki, Stirling 



Heavy sfermions 



• In the first step we use only the masses of xf, Xl^ 
X2 and the chargino pair production cross section, in- 
cluding the full leptonic and the semileptonic decays 
as observables. A four-parameter fit for the parame- 
ters Ml, M2, fi and mj, has been applied. 

• In the second step we include as an additional observ- 
able the leptonic forward-backward asymmetry. Only 
the semileptonic and purely leptonic decays were used. 
The SU(2) relation between the two virtual masses mc, 
and rrie^ has been applied as an external constraint. 

• As an attempt to test the SU{2) mass relation for the 
slepton and sneutrino masses, in the last step both the 
leptonic and hadronic forward-backward asymmetries 
have been used. A six-parameter fit for the parameters 
Ml, M2, /i, rriiy, rriei^ and tan/3 has been applied. 

a) Analysis without forward backward asym- 
metry 

We use as observables the masses m-±, m^o and 
the polarized chargino cross section multiplied by the 
branching ratios of semileptonic chargino decays, see 
Section O and Tabled 

We apply a four-parameter fit for the parameters 
Ml, M2, fi and m^^ for fixed values of tan/? = 5, 10, 
15, 20, 25, 30, 50 and 100. Fixing tan/3 is necessary for 
a proper convergence of the fitting procedure because 
of the strong correlation among parameters [10]. We 
perform a test. It turns out that for tan/3 < 1.7 the 
measurements are inconsistent with theoretical predic- 
tions at least at the la level. We obtain: 

59.4 < Ml < 62.2 GeV, 
118.7 < M2 < 127.5 GeV, 

450 < fi < 750 GeV, 
1800 < m£>^ < 2210 GeV. 

Figure [1] shows the migration of la contours in rrij^^- 
M2 (left), M2-IJ. (middle) and M1-M2 (right), the other 
two parameters being fixed at the values determined 
by the minimum of for tan/3 changing from 5 to 10, 
20 and 50. 

b) Analysis including leptonic forward 
backward asymmetry 

We now extend the fit by using as additional observ- 
able the leptonic forward-backward asymmetry for po- 
larized beams (-90%, -1-60%). 

As a result the multiparameter fit strongly improves 
the results. No assumption on tan/3 has to be made in 
the fit since for too small or too large values of tan/3 
the wrong value of ApB is predicted. We find 

59.7 < Ml < 60.35 GeV, 
119.9 < M2 < 122.0 GeV, 
500 < fi < 610 GeV, 
1900 < m£>^ < 2100 GeV, 
14 < tan/3 < 31. 

The constraints for the mass nio^ are improved by a 
factor of about 2 and for gaugino mass parameters Mi 



and M2 by a factor of about 5, as compared to the 
results of previous section with unconstrained tan/3. 
The error for the higgsino mass parameter /i also de- 
creases significantly. Taking the constraints of eq. pOI) 
leads to a prediction of 

506 < m^o < 615GeV, 
512 < m^o < 619GeV, 
514 < TO-± < 621GeV. 

c) Analysis including hadronic and leptonic forward- 
backward asymmetries: test of SU{2) 

With the constraints for the squark masses from the 
LHC, the hadronic forward-backward asymmetry could 
be used to control the sneutrino mass. The leptonic 
forward-backward asymmetry provides constraints on 
the selectron mass and the SU(2) relation could be 
tested. 

We perform a scan of the parameter space, taking 
light chargino and neutralino masses, 4 cross section 
and leptonic asymmetry measurements and apply a 
test. We derive the following constraints: 

59.30 < Ml < 60.80 GeV, 
117.8 < M2 < 124.2 GeV, 
420 < fi < 950 GeV, 
1860 < m£>^ < 2200 GeV, 
1400 GeV < mg^, 
11 < tan/3. 

Including hadronic forward-backward asymmetry im- 
proves the constraints as follows: 

59.45 < Ml < 60.80 GeV, 
118.6 < M2 < 124.2 GeV, 
420 < fi < 770 GeV, 
1900 < m£>^ < 2120 GeV, 
1500 GeV < mg^, 

11 < tan/3 < 60. 

The most significant change is for the sneutrino mass, 
for which error bars become smaller by ~ 50%. Also 
an upper limit on tan/3 is found. However we do not 
get an upper limit for the selectron mass. Nevertheless, 
the results for the selectron and sneutrino masses are 
consistent with the SU{2) relation. 



4 Conclusions 

Scenarios with heavy scalar particles seem to be very 
challenging for determining the MSSM parameters since 
only a very limited amount of experimental informa- 
tion can be accessible. 

A very powerful tool in this kind of analysis turns 
out to be the forward-backward asymmetry. This asym- 
metry is strongly dependent on the mass of the ex- 
changed heavy particle. If the SU{2) constraint is ap- 
plied, the slepton masses can be determined to a pre- 
cision of about 5% for masses around 2 TeV at the 
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Table 1. Masses of the SUSY particles (in GeV). 
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Fig. 1. Migration of la contours for tan/3 — 5, 10, 20, 50 
with the other two parameters fixed at the values deter- 
mined by the minimum of for each tan (3 |10] . 



Table 2. Branching ratios for some important decay modes 
in the studied MSSM scenario, £ = e,/i, r, g„ = u,c, Qd = 
d, s. Numbers are given for each family separately. 



Mode 


9 ^ X2bb 


5 -* Xi Quqd 


Xt XildQu 


BR 


14.4% 


10.8% 


33.5% 


Mode 




ii,2 -* gt 


Xl-^xlrvt 


BR 


3.0% 


66% 


11.0% 



Table 3. Cross sections for the process e^e~ Xi XT 
and forward-backward asymmetries (Afb) in the leptonic 
X\ Xi^~'^ s-iid hadronic Xi ~^ Xi'SC decay modes, for 
different beam polarization P^- , P^+ configurations at cms 
= 350 GeV and 500 GeV at the ILC. Concerning the 
errors, see text. 







{P.-,P.+ ) 


75 = 350GeV 


(-90%, +60%) 


(+90%, -60%) 


'^ixtxT)/^ 


6195.5 


85.0 




1062.5±4.0 


14.6+0.7 


Afb (/-)/% 


4.42±0.29 




Apb(c)/% 


4.18±0.74 




= 500GeV 


(-90%, +60%) 


(+90%, -60%) 


f^(x^xr)/fb 


3041.5 


40.3 


l^(X^Xr) Bslc Gaic/fb 


521.6+2.3 


6.9+0.4 


AFB(r)/% 


4.62+0.41 




>1fb(c)/% 


4.48+1.05 





ILC running at 500 GeV. Although only very limited 
information is available, powerful predictions for the 9 
heavier charginos / neutralinos can be made. 
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